Description of rank four PPT entangled states of two qutrits 
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It is known tliat some two qutrit entangled states of rank 4 with positive partial transpose (PPT) 
can be built from the unextendible product bases (UPB) [ Phys. Rev. Lett. 82, 5385 (1999) ]. We 
show that this fact is indeed universal, namely all such states can be constructed from UPB as 
conjectured recently by Leinaas, Myrheim and SoUid. We also classify the 5-dimensional subspaces 
of two qutrits which contain only finitely many product states (up to scalar multiple), and in 
particular those spanned by a UPB. 
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I. INTRODUCTION 

The positive-partial-transpose entangled states (PPTES) arc of particular importance and interest in quantum 
information (for a review see [11]). For a state p acting on the Hilbert space % :— Ha ^Hb, the partial transpose 

computed in an orthonormal (o.n.) basis {|ai)} of system A, is defined by :— We say that p 
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is a PPT [NPT] state if p'" > [p"" ^ 0, i.e., has at least one negative eigenvalue]. The most intriguing feature 
of PPTES is their non-distillability under local operations and classical communications (LOCC) [IJ]. This means 
that a PPTES, say p, cannot be locally converted (asymptotically) into a pure entangled state even if infinitely 
many copies of p are provided. Since most quantum-information tasks require pure entangled states, a PPTES is 
a quantum resource which cannot be directly used in practice [1]. Nevertheless in the past few years the PPTES 
have been extensively studied in connection with the phenomena of entanglement activation and universal usefulness 



[25l . l33l | , the distillable key [22| , the symmetry permutations [38[ and entanglement witnesses [32| , both in theory and 
experiment. 

For many applications, it is an important and basic problem to decide whether a given PPT state is entangled 
or separable, i.e., the convex sum of product states \ai,hi){ai,bi\ (39j . The separability problem has an extensive 
application in quantum information, metrology, computing, quantum non-locality, and mathematics (like positive 
maps and C*-algebras). Moreover, the problem has been proved to be NP-hard and hence it attracted a lot of 
attention from computer scientists [13]. In 1996, the first necessary condition was given by Peres [13], saying that 
the separable states are always PPT. So to solve the problem it suffices to consider only the PPT states. Next in 
1997, the Horodeckis [1^ showed that this is necessary and sufficient for 2 (g) 2 and 2 (g) 3 states. However both of 
these cases lack the PPTES. Actually since the first PPTES was constructed [13], researchers lacked for a long time 
the analytical characterization of PPTES in any bipartite systems of given rank and local dimensions. For example, 
it was surprisingly difficult to decide whether a given state of rank 4 in 3 (g) 3 space is a PPTES, which is also the 
smallest space in which PPTES may exist [1^. A well-known method [^ for construction of PPTES proposed in 1999 
was based on the unextendible product bases (UPB). It is applicable to arbitrary bipartite and multipartite quantum 
systems. Another systematic method for two-qutrit systems was provided by Chen and Dokovic who proved in 2010 
that a PPT state of rank 4 is entangled if and only if there is no product state in its range [3 • 

Our main result (see Theorem I25p shows that any two-qutrit PPTES p of rank 4 can be constructed from an 
unextendible product basis (UPB) [J by using the method proposed by Leinaas, Myrheim and Sollid in [s^- We state 
their conjecture formally as Conjecture [TS] Explicitly, we prove that (up to normalization) p = A® B n{-!/)} ® B\ 
where {?/'} is a UPB, n{-0} is the projector on the subspace orthogonal to {■(/'}, and A,Be GL3, the group of invertible 
complex matrices of order 3. Let us point out that the papers [20l . [sj provided strong numerical evidence for the 
validity of this result and motivated us to pursue this study. Moreover, the authors of these papers suggest that the 
higher dimensional PPTES may have similar properties, although some of them have to be modified. This may be of 
interest for further research in this direction. 

It is well-known that the set iSsep of separable states and the set iSppx of PPT states are both compact and convex, 
and that 5sep ^ ^ppx. A basic task is then to characterize their extreme points. The latter are the states which are 
not convex sums of other states in the convex set. It is also a well-known fact that we can generate any state in such 
a set by taking the convex sum of extreme points. Though it is known that the extreme points of 5sep are exactly 
the pure product states [3§| , we know quite little about the extreme points of iSppt [l^l • Only in the case of 2 (g 4 
systems, the PPTES have been partially classified by the ranks of both the state and its partial transpose [l[. Here, 
we will show that all two-qutrit PPTES of rank 4 are such extreme points (and also edge PPTES). 

We also show that no PPTES of rank 4 exist in the symmetric subspace of the two-qutrit system. Furthermore, it 
is known that the UPB basis states are not distinguishable under LOCC [2|. Thus we exhibit the essential connection 
between the LOCC-indistinguishability and PPTES. 

The content of our paper is as follows. In Sec. |TT]we define the notion of general position for m-tuples of product 
vectors (or points on the Segre variety x V^), see Definitions [T] and HI We also define and study the properties 
of the biprojective (BP) equivalence of such m-tuples, and in particular for quintuples. To a quintuple of product 
vectors in general position we assign six GL3 x GL3 invariants {J.f', JP), (« = 1, 2, 3), and show that two quintuples are 
BP-equivalent if and only if they share the same invariants. These invariants are subject to the relations J^^J^J^ = 1 
and JiJiJi = 1, and can take any complex values except and 1. We then investigate the product states contained 
in the 5-dimensional subspace spanned by a quintuple of product states in general position. In particular, it follows 
form the part (a) of Proposition [7] that if p is a two-qutrit PPT state of rank 4 then its kernel contains only finitely 
many product vectors (up to scalar multiple). This fact is used later to prove Theorem 1221 

In Sec. mil we consider the 8-dimensional complex projective space associated to the 3 (g 3 Hilbert space T-L. We 
consider the case where a 4-dimensional projective subspace, T'^, and the Segre variety £2,2 intersect only at finitely 
many points. Then there are at most 6 points of intersection and to each of them one assigns (in Algebraic Geometry) 
a positive integer known as the intersection multiplicity. The sum of these integers is necessarily equal to 6. We refer 
to the corresponding partition of 6 as the intersection pattern of 7^"* and £2.2- Wc show that all feasible patterns, i.e., 
all partitions of 6, occur. For instance, there exists a 5-dimensional vector subspace of % which contains only one 
product state (up to a scalar multiple). 

In Sec. IIVI we study the invariants of quintuples of product states in general position which arise from UPB. These 
invariants are real and take values in one of the open intervals N (— oo,0), p (0,1) and P (l,+oo). By 
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replacing the values of the 6 invariants by the letter designating the interval to which the invariant belongs, we obtain 
the 6-letter symbol attached to the quintuple. We show that among the 120 permutations of the 5 product states of 
a UPB, there are exactly 12 different symbols that arise in this way. We say that these 12 symbols are UPB symbols. 
The same assertion is valid if we use all 6 product states contained in the 5-dimensional subspace spanned by a UPB 
and construct from them all 720 possible quintuples. Their symbols are also UPB symbols. This is the key tool in 
the proof of our main result. 

In Sec. |V] we prove that the kernel of any PPTES p of rank 4 contains exactly 6 product states (up to scalar 
multiple) and that these 6 states are in general position. Next we show that there exist A, B S GL3 such that the 
transformed state a :— A® B p ® B'^ \s invariant under partial transpose, i.e., = a. We also show that if two 
normalized PPTES of rank 4 have the same range (or kernel) then they are equal. Finally, by making use of the 
UPB symbols we prove our main result. Theorem [55] We conclude that section with a description of the stabilizer Gp 
of a PPTES p of rank 4 in the group PGL := PGL3 x PGL3. We show that this is a finite group isomorphic to a 
subgroup of the symmetric group Sq of degree 6, and we compute these groups for the two most interesting and most 
symmetric cases arising from Pyramid and Tiles UPB. 

In See. IVII we discuss a few physical applications of our results. We show how to build analytically the PPTES of 
rank 4 of two qutrits. This is also helpful for the construction of entanglement witnesses which detect these PPTES. 

Throughout the paper we use the expression "A/ x N state" to mean a bipartite state, say p, whose local ranks are 
rank pA = M and rank ps = N. By TZ{p) we denote the range of a linear operator p. 

Another proof of the Conjecture [TSl was obtained by L. Skowronek |35(. 

II. PRODUCT STATES OF TWO QUTRITS 

For convenience, we shall represent any product vector \x,y) as a 3 x 3 matrix [xii/j] where i,j S {1,2,3}. In 
this notation, product vectors correspond to matrices of rank 1. We shall often work with these vectors only up to 
scalar multiple, in which case we consider them as points in the 8-dimensional complex projective space associated 
to the 9-dimensional Hilbert space H = Ha ® Hb- Then the product vectors form the so called Segre variety 
^2,2 y<V^ Q V^, isomorphic to the product of two copies of the complex projective plane V^. 

A. Projective invariants Ji, J2, J3 

The complex general linear group in dimension 3, that is the group of invertible complex matrices of order 3, will 
be denoted by GL3. We shall also use the group GL := GL3 x GL3, the direct product of two copies of GL3. It acts 
naturally on H = Ha ®T~Lb via invertible linear operators (ILO) A® B. The points of the complex projective plane 
will be identified with the 1-dimensional subspaces of Ha (or He)- As GL3 permutes the 1-dimensional subspaces 
of Ha, it induces an action on . The subgroup C^, consisting of nonzero scalar matrices of GL3, acts trivially on 
7-"^. Thus we obtain an action of the quotient group PGL3 := GL3/C^ on . This group is known as the complex 
projective general linear group (in dimension 3) and the transformations that it induces on are known as projective 
transformations. The direct product PGL PGL3 x PGL3 acts on the Segre variety x , and we refer to it as 
the group of biprojective transformations of S2,2- 

Invertible linear operators are also useful in quantum information. Two n-partite states p and a are equivalent under 
stochastic LOCC (or SLOCC- equivalent) if p = {g)"^^ A,cr{g)"^i A\ for some ILO Ai, . . . , An. They are LU- equivalent 
if Ai can be chosen to be unitary. In most cases of the present work, we will have n = 2. Both LOCC and SLOCC are 
referred to as physical operations in quantum information [23. The essential difference between them is that LOCC 
can be implemented with certainty while SLOCC succeeds only with some nonzero probability. 

Definition 1 We say that an m-tuple of pure states (I'/'fc)))!^^ ('^^ Hb) is in general position if any two or 

three of them are linearly independent. If {\4>k))^^Q and {\4''kl)l^=o "^^^ ^'^'^ m-tuples in Ha such that A\(l>k) oc 
for some invertible matrix A, then we say that they are projeetively equivalent or P-equivalent. 

Let us recall the following elementary fact from Linear Algebra also known as the Four Point Lemma (see e.g. [T^ . 
Lemma 11.2]): 

Lemma 2 // (|(?ifc))Lo '^^id (|0fe))| are quadruples of pure states in general position in Ha, then there exists an 
invertible matrix A, unique up to a scalar factor, such that A\(j)k) oc for each k. 
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In particular, for any quadruple {\4'k))k=o general positfon there exists an invertible matrix A such that ^|^fc) cx 
|fc) for A: = 0, 1, 2 and A|(/)3) = |0) + |1) + |2). For convenience, we say that A transforms the quadruple {\4'k))k=o '^^^^^ 
the canonical form. One can construct A explicitly as follows. Let X — [\(j)o) |02)] and let D be the diagonal 
matrix whose diagonal entries are the components of the vector X~^\(j>s). Then we have A = D~^X~^. 

In the language of Projective Geometry the above lemma can be restated as follows. If {Pk)k=o and (^fc)fc=o 
two quadruples of points in the complex projective plane in general position (i.e., they are distinct and no three 
points Pfc are colinear, and similarly for the points P^) then there is a projective transformation T of such that 
r(Pfc) = PI for each k. 

However, the action of GL3 on quintuples of points in general position in V'^ is not transitive. Indeed let (P^ = 
I0fe))fc=o such a quintuple. Then all determinants 

Aij,fe = det[|(?!)i) |0fc>], distinct), (1) 

are nonzero. The rational functions 

J ^2, 0,4^0. 1,3 J _ ^0. 1.4^1.2,3 Ai_2.4A2,0.3 /r,\ 

^ ~ 7^^ A ' 2 - TT T , J3 - -7 -T i^j 

^2,0,3^0,1,4 ^0,1,3^1,2,4 ^1,2,3^2,0,4 

are projective GL3-invariants of quintuples in general position. These invariants may take arbitrary complex values, 
except and 1, subject to the relation Ji J2J3 = 1. The following result follows easily from the Four Point Lemma 
and the fact that P4 is determined uniquely by the quadruple (Pfc)|^g and the values of the invariants Ji. 

Proposition 3 Two quintuples of points (P^ = \4>k))\^o o.nd (P^ = \'t''k))t=o general position are P-equivalent if 
and only if they share the same values of the three invariants Ji . 

This means that if the two quintuples satisfy the invariance conditions, then there exists A g GL3 such that 
A\(t>k) — Cfe|0fe) for some scalars c^, which may be all different. 

Definition 4 We say that an m-tuple (iV'fe) \4'k) ® \Xk})^Zo non-normalized product states m a 3 x 3 system is 
in general position if each of the m-tuples {\(t>k))k^J^o ^'^'^ (IXfe))fc^^ ^■s ''^ general position. In that case we also say 
that the corresponding m-tuple of points (Pfc = |'0fe))fc!=o^ on 5]2,2 is in general position. We also say that two m-tuples 
of product states (|V'fe) = \(pk) ® \Xk))T=o "-^^ (|Vfc) = Wk) ® \x'k))T= 

Q are biprojectively equivalent or BP-equivalent 

if there exists A(E} B G GL such that {A (E) B)\tpk) oc 1-0^) for each k. 

We shall use the same terminology for the m-tuples of points lying on £2,2- 

In the important case to = 5, we have two sets of invariants, one for {\4>k))t=o ^^'^ ^he other for (|xfe))t=o- 
shall denote the former by and the latter by jf . 

The following proposition is an immediate consequence of the definitions and results stated above. 

Proposition 5 The group GL acts transitively on the quadruples of points in general position in ^2,2. Two quintuples 
of points {Pk = |</)fc) ® |Xfc))Lo a"*^ = Wk) ® lx'fe))i 

in Tj2.i, both in general position, are BP-equivalent if and 
only if the quintuples of points {\(j)k))'k=o ^''^d i\4>'k))k=o in-P^ '^'"6 P-equivalent and the same is true for the quintuples 
of points (|xfc))Lo and i\x'k) )t^o. 

We apply our results to the manipulation of families of quantum states. That is, we consider the condition on which 
two sets of quantum states {pi, . . . , p„} and {tri, . . . , (T„} can be (probabilistically) simultaneously convertible via a 
physical operation e, i.e., e{pi) :— Aj{pi)A}j = Oi, for all i Q. The problem is in general difficult even for the case 
of single-party and n = 2, which has been studied in terms of single-party states for distinguishing quantum operations 
(1^ . In present work, pi and Ui are single-party pure states in general position. To realize the conversion, the Kraus 
operators Ai have to be pairwise proportional. So the Four Point Lemma and Proposition [3] can be used to decide 
the simultaneous conversion between two sets of 4 and 5 qutrit states in general position, respectively. Moreover, 
Proposition [5] works for the conversion between two sets of bipartite product states. Hence we have produced some 
new manipulatable families of states. Further, we have 

Lemma 6 // (|'0A;))|^g is a quadruple of product states in general position, then the subspace that they span contains 
no other product state. 
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Proof. Let l^/jfe) = \(f)k) ® \Xk), fc = 0, . . . , 3. By applying the Four Point Lemma in T-La and T-Lbi we may assume 
that the quadruples (|'0fc))fe=o canonical form 



iVfc) (X |/c/c), (fc = 0,l,2), IV'a) (X 



1 1 1 
111 
111 



(3) 



Note that the space of 3 x 3 diagonal matrices contains only 3 product vectors (up to a scalar multiple). Hence the 
assertion of the lemma follows from the fact that a non-diagonal matrix 



a 5 5 
6/3 6 
S S J 



has rank 1 if and only iia~(3 = "f = S. 



(4) 



□ 



B. Intersection of and E2,; 



Projective geometry has quite different properties from those of the affine geometry. For instance, in the complex 
afhne plane, C^, there exist distinct straight lines which are parallel (and so do not meet). However, in the complex 
projective plane, P^, any two distinct lines meet at exactly one point. More generally, any two projective varieties X 
and Y in 7"* must meet if DimX + DimF > n. 

Let V he a 5-dimensional subspacc of the 3x3 system H = T-La®T~Lb and the projective 4-dimensional subspace 
associated to V . Since the complex dimensions of and ^2,2, namely 4 and 4, add up to the dimension 8 of the 
ambient projective space , these two varieties must have nonempty intersection, i.e., V always contains at least one 
product state, see e.g. f2ll . Proposition 11.4]. On the other hand V may contain infinitely many product states, i.e., 
the intersection of this V'^ and the Segre variety ^2,2, may have positive dimension. Let us assume that V contains 
only finitely many product states which we treat as points in , say Pi, i = 1, . . . This is often expressed by 
saying that in this case the varieties V'^ and S2,2 intersect properly. Since and I]2,2 have degrees 1 and 6 (see [2l|, 
p. 233]), respectively, Bezout Theorem tells us that 1 < < 6. More precisely, we have /^i + • • • + /i^ = 6 where 
/ii, a positive integer, is the intersection multiplicity of the point Pi. See e.g. |l3. [2l| for more details about these 
multiplicities and the Bezout Theorem. 

We next consider an arbitrary quintuple {\ipk))t=o product vectors in general position and the 5-dimensional 
subspace V that they span. We would like to determine whether V contains any additional product states. The 
following proposition gives the answer. Since the invariants determine uniquely (up to BP-equivalence) the quintuple 
of points on £2,2 hi general position, it is possible to analyze and answer the above question solely in terms of the 
invariants. 



Proposition 7 Let (Pfe = {ipk) = |<^fe) O \xk))t= 

J,P [i — 1,2,3) he its invariants. Denote by V the subspace spanned by the jV'fc 
space. We consider the five equations: 



be a quintuple of product states in general position, and let jf and 

and by the associated projective 



(i = 1,2,3), 



J2^(l - Jt){l - J^) = J^{1 Jf )(1 - J2^), 
J^l - J2^)(l - Jf ) = Jf (1 - J2^)(l - J^). 



(5) 
(6) 

(7) 



(a) If any two of these equations hold, then all of them do. In that case V contains infinitely many product states. 
Moreover, any state p with ker p = V must be NPT. 

(b) If exactly one of the above equations holds, then V contains no additional product states. 

(c) If none of the above five equations holds, then V contains exactly one additional product state, say \^) = |c, z) = 
[ciZj]. The vectors \c) and \z) are given (up to a scalar factor) by the formulae 



(1 - JF)/iJf - ^1^) 

J2^(l- J2^)/(J2^- J2^) 
(1- jf)M^(jf - J^) 



(8) 



Ji{l-Ji)/{Ji-Ji) 
{l-Ji)/J?{Ji~Ji) 



(9) 
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where A and B are the matrices which transform the quadruples (|<^fe))fc=o and (|xfc))fc=o canonical form, 

respectively. Moreover the six product states \tpk), = 0, . . . ,4 and {tp) are in general position. 

Proof. We begin by transforming the quintuples {\4'k))t^Q and (|xfc))fc=o matrices A and B, respectively. Thus 
we may assume that the quadruples (|'?I>fc))^^o ^^"^ (IXfe))fc=o the canonical form, i.e., 



|V'fc)«|fcfc> (fc = 0,l,2), 
IV'4) = \b,y) [hyj]. 



(10) 
(11) 
(12) 



Since the Pk are in general position, all components hi and pi are nonzero, and hi ^ bj and yi ^ yj for i ^ j. A 
computation shows that 



J? 



^2/63, 4^b^lbu J^^b,/b2; 

2/2/2/3, J2 = 2/3/yi, J3 = yi/y2- 



Thus we have ^ 1 and Jf 7^ 1 for z = 1, 2, 3, and the equations © and ([7]) can be written as 

biy2 + b2y3 + bsyi = 612/3 + ^2^1 + 63y2, 
111 111 



6i2/2 &22/3 ^32/1 ^12/3 ^22/1 ^32/2' 



(13) 
(14) 



(15) 
(16) 



respectively. 

We consider first the case (a). If two of the equations (O hold, so does the third because of the identity J1J2J3 = 1- 
Clearly, the remaining two equations also hold. The other cases can be treated similarly. Let us just consider the 
hardest case where the two equations displayed above hold. By solving the first for y^ and substituting into the 
second, we obtain that 



(2/1 - 2/2)(fei - &3)(^2 - b3){biy2 ~ 622/1) 
bib2b3yiy2{biy2 - ^22/1 + bayi - 632/2) 



(17) 



Note that 612/2 — 622/1 + 632/1 — 632/2 7^ because 2/3 7^ 0. Hence we conclude that 612/2 — 622/1 = 0, i.e., J3 = J3 . This 
means that this case reduces to one of the other cases, and the first assertion of (a) is proved. 

We now assume that all five equations hold. We can further assume that hi = yi = 1 and, consequently, 62 = 2/2 
and 63 = 1/3. Thus V consists of all symmetric matrices 



u a + /362 a + /363 
a + I3b2 V a + j3h2hz 
a + /363 a + /36263 w 



u, V, w,a, P ^ C. 



By specializing the diagonal entries 

(a + /362)(a + /363) 

u = . 

a + P6263 



(a + ^62)(a + /362 63) 
a + /363 



{a + Pb3){a + /3b2b3) 
a + Pb2 



(18) 



(19) 



in the above matrix (|18p , we obtain a family of non- normalized product states depending on two complex parameters 
a and j3. Since it is contained in V , the second assertion of (a) is proved. 

Let p be a non-normalized state with kerp = V. Its range, V^, is spanned by |01) — |10), |12) — |21), |20) — |02) 
and \ip) := 6(|01) + |10)) - (1 + 6)(|12) + |21)) + (|20) + |02)), where b* = -63(1 - 63)762(1 - 63). Without loss of 
generality, wc can write p — X]?=o where \6i) are linearly independent non- normalized pure states 

%) = |01)-|10), (20) 

= ao(|01)-|10)) + ai(|12)-|21)), (21) 

1^2) = a2(|01) - |10)) + a3(|12) - |21)) + a4(|20) - |02)), (22) 

1^3) = 05(101) - |10)) + a6(|12) - |21)) + a7(|20) - |02)) + x\^), (23) 
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with ai,a4,x > 0. Assume that (t := > 0. Since |00) G kerp, the first diagonal entry of the matrix a is 0. 
Consequently, its first row must vanish. Hence we obtain that 



al + {x — a7){x + a^) = 0, 

a2a4 + {bx + a5){x + aj)* = 0, 

0304 + (x — a7)(fe.T — 05)* = 0, 

-1 - \ao\^ - \a2\^ + {bx + a5){bx - as)* = 0. 



we deduce that 07 is real and 



From Eq. 

02 = 6(07 — a;^)/a4, and then Eq. (|27p gives the contradiction —1 — |aol'^ = 



From the next two we deduce that 05 
2 



(24) 
(25) 
(26) 
(27) 

-baj and 



Hence our assumption that a > must be false, i.e., p must be NPT. Thus all three assertions of (a) are proved. 
Next we consider the case (b). Assume that one of the Eqs. ([S]) holds. Say, = , i.e., 61^2 — &2yi = 0. Then 
we can also assume that &i = yi = 1, and so 62 = 2/2- Hence V consists of all matrices 



X 



u a 

f /3fo2 



f /3&2 
V 

Pb2b3 



- ^^22/3 

w 



u, V, a, /3 e C. 



(28) 



Assume that such a matrix X has rank 1. Let us also assume that a/3 7^ 0. Then we can assume that (3 = 1. The 
following equations must hold 



uv = {a + 62)^, 

u{a + 6263) = (a + b2)ia + 63), 
w(a + 622/3) = (a + 62)(a + Va)- 



(29) 
(30) 
(31) 



From the last two equations we obtain that u{a + b2b3){a + 2/3) = + 622/3)(a + 63), i.e., a(&2 — 1)(^3 ~ 1/3)^ = 0- As 
we deal with case (b) and we assumed that = and 61 = ?/i = 1, we must have 7^ J^, i.e., 63 7^ 2/3- Hence, 
we deduce that u = 0, and so a = —62- Since X has rank 1 and 62 7^ 63, we deduce that the (2,3) entry of X must 
vanish. This gives the contradiction 62(2/3 — 1) = 0. Consequently, we must have a/3 = 0. Then Lemma [S] implies that 
X oc lipk) for some fc = 0, . . . , 4. 

The case when Eq. (|6]) holds can be reduced to the above case by switching the states IV'2) and \ip3). After this 
transposition, the new values of the invariants J^^ are given by the formulae: 



J[=l- 



j,A_ 
J2 — 



•-'2 



•^2 



1' 



J3 — 



J2^-l 
J2^(l- Jl^)' 



(32) 



J3 is the same as ©■ 



and the same formulae are valid for . It remains to observe that the equality J3 

The case when Eq. ([7]) holds can be reduced to the previous case by switching the states IV'3) and \ip4)- After this 
transposition, the new values of the invariants J^^ are just the reciprocals of the old invariants J/^, and similarly for 
J'^^ . Then Eq. ^ is replaced by Eq. ([7]). This completes the proof of (b). 

It remains to prove (c). This is a generic case; the five subcases of (b) can be intuitively viewed as the limiting 
cases of (c) which occur when the state (which we are going to construct) becomes equal to one of the given five 
states so that its multiplicity increases from 1 to 2. 

Since \ih) G V, we have 



(33) 



Lemma ini implies that all coefficients must be nonzero. The six off-diagonal entries of these matrices give the following 
system of equations 



C1Z2 = a + /36i2/2, C1Z3 = a + /36iy3, C2Z3 = " + /3622/3, 
C2Zi = a + /3622/1, C3Z1 = a + /363yi, C3Z2 = « + /363?/2- 



(34) 
(35) 



By eliminating ci from the first two equations, and C2 and C3 from the last four, we obtain a system of linear equations 
in the unknowns z,-: 



(a + /36iy3)z2 - (a + ^612/2)2:3 = 0, 
-(a + /3622/3)2i + (a + /3622/i)^3 = 0, 
(a + /363y2)zi - (a + /3632/i)2:2 = 0. 



(36) 
(37) 
(38) 
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Since this system has a nontrivial solution, its determinant must vanish: 



a + l3biy3 -(a + /3&iy2) 

-(a + /3&2J/3) a + l3b2yi 

a + fib^yi -(a + /3632/1) 



= 0. 



This gives the equation 

(a + I3biy2){a + I3b2yz){a + Pb^yi) = (a + I3biys.){a + I3b2yi){a + /3&32/2)- 

After expanding both sides, the terms with and f3'^ cancel and after dividing both sides by a/3, we find that 

hhyiiva - 2/2) + &i632/2(yi - Vs) + &i&2?/3(y2 - Vi) 
biiys - V2) + &2(yi - Vs.) + b'i{y2 - vi) 



a 






cx 



(39) 



(40) 



(41) 



We can replace the proportionality sign with the equality because |c) and \z) are determined only up to a scalar factor. 
The system of linear equations p6p for the Zi has the unique solution up to a scalar factor: 



\z) (X 



(a + I3b2yi)ia + pb^yi) 
{a + I3b2yi){a + /363J/2) 
(a + /362y3)(a + pbsVi) 



(42) 



where a and /3 are given by Eq. (j41[) . After substituting the expressions for a and P and cancelling two factors, we 
obtain the formulae From ([55)1 we have 



ziZ3(a + /361J/2) 
ZiZ2(Q! + I3b2y3) 
Z2Z'i[a + /3&3yi) 



(43) 



and, by using (PT|) . we obtain ([5]). 

Finally by using the above expressions, we can verify that the six product states IV'fc), — 0, • • • ,4 and j'i/') are 
indeed in general position. This completes the proof. □ 

We remark that, by Bezout Theorem, in the case (b) exactly one of the 5 intersection points of n S2.2 must have 
multiplicity 2 and the other multiplicity 1. If jf" = Jf holds, then the point with multiplicity 2 is Pi-\. If Eq. © 
[Eq. dll)] holds, then this is the point P3 [P4]. 

The case (c) is of special interest and we single it out in the following definition. 

Definition 8 A quintuple of product states {\'4'i))j=o is regular if it is in general position and the 5- dimensional 
subspace spanned by the \ipi) contains exactly one additional product state (up to scalar multiple). 

As an immediate consequence of the above proposition, we observe that in the case (c) the map sending \y)) — >■ 
(|c), \z)) is involutory, i.e., it also sends (|c), \z)) — >■ (|&), \y)). As another consequence, we have 

Corollary 9 Let (|'0i))i=o be a regular quintuple of product states and its invariants. Denote by lip) the 

unique additional product state in the subspace spanned by the \ipi). Then the invariants J[ , J,' of the quintuple 
(iV'o), iV'i), l'02), |'/'3), 1^)) are given by the formulae 



{i-Ji){Ji-Ji) 



Ji{i~Ji){J^ 



Jim 



jH^~j?){Ji-Jiy 



•''3 



•''3 



{i-j?){Ji-Ji) 

Ji{i-Ji){jF-J^)' 



(44) 
(45) 



All cases (a-c) of the above proposition may occur; for (b) and (c) sec the last two examples in the proof of Theorem 
1101 In particular, the proposition shows that the number of product states in V may be infinite or only 5 even when 
we impose the condition that the 5 given product states are in general position. Nevertheless, we will show later that 
the kernel of a 3 x 3 PPTES of rank 4 always contains 6 product states, and that they are in general position (see 
Theorem [221 in Sec. |V|. Hence, if the kernel of a state p of rank 4 is of type (a) or (b) of Proposition [3 then p must 
be NPT. This may shed new light on the problem of entanglement distillation of 3 x 3 NPT states of rank 4 [7|. We 
investigate further the properties of 5-dimcnsional subspaccs in the next section. 
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The results proved in this section may help to solve another long-standing quantum-information problem, namely 
the state transformation under SLOCC [3] . The problem is completely solved for 2 x M x iV pure states 0, § ; however 
it becomes exceedingly difficult when all three dimensions are bigger than two, e.g., deciding the SLOCC-equivalence 
of two 3x3x3 states. Here we consider the transformation between two 3x3x5 states G 'Ha®T~Lb®T~Lc and 

assume that there are 5 product states \ai,hi) £ Tyq in general position. According to the range criterion @, 

IV'), If) are interconvertible via ILOs Va, Vb, Vc if and only if Va ® Vsjai, hi) G TZ{(pab) where ipAB = Trc |<y3)(¥'|. By 
finding out the product states in TZ{(pab), wc can decide the SLOCC-equivalcncc of \(p) via Proposition [5] and [T] 
Hence, the SLOCC-equivalence of two 3x3x5 states can be operationally decided provided that 5 product states in 
their range of bipartite reduced states are available. Furthermore we can expand subspaces, such as the 3x3x3 and 
3x3x4 subspaces by respectively adding 2 or 1 linearly independent (product) states, to span the whole 3x3x5 
space. Thus we may treat the SLOCC-equivalence of tripartite states of the former subspaces similar to the latter, 
when we can build 5 product states in the range of corresponding reduced states. 



III. PRODUCT STATES IN 5-DIMENSIONAL SUBSPACES 



As in the previous section, let us consider the intersection V* H ^2.2 and assume that it is proper. Our main 
objective here is to investigate various possibilities for this intersection and provide concrete examples for each case. 



A. Intersection patterns 

Recall that because Y.2,2 is a finite set, consisting say of k points, we know that necessarily A: < 6. Denote by 
fii the intersection multiplicity of the point Pi . When arranged in decreasing order /ii > ^2 ^ • • • > Mfe , they form a 
partition (^1, . . . , /i^) of the integer 6. We shall refer to this partition as the intersection pattern (see [l4l p. 182]). 
Altogether there are 11 such partitions and we shall first prove that all of them occur as intersection patterns. 

Theorem 10 All 11 partitions of 6 occur as intersection patterns ofV^ H ^2,2, where is a complex ^-dimensional 
projective space and T.2.2 the Segre variety. 

Proof. We just have to provide examples of two qutrit states p of rank 4 whose kernel, viewed as a complex projective 
4-dimensional subspace, has the specified partition of 6 as its intersection pattern with T,2,2- We shall subdivide the 
list into 6 parts corresponding to the number, say k, of product states contained in ker p. The corresponding partitions 
of 6 are those having exactly k parts. For k = 6 the examples will be provided by the normalized projectors associated 
with the UPB (see Theorem and also by generic separable states of rank 4 (see Lemma ETjl . However we shall 
include a concrete example with fc = 6 in our list. 

In each case we assume that p = X]i=o l^j)(^j| and give the formulae for the pure states \ipi). We also list the k 
product states in the ker p as well as their intersection multiplicities pi. These multiplicities were computed by means 
of the free software package Singular [l5| for symbolic computation in Commutative Algebra. 

For fc = 1 we set 

|Vo) = |12), |V^i)-|21), |V^2) = |01)-|10)-|22), IV-s) = |02) + |11)-|20). (46) 

The kernel is spanned by the states |00), |01) + |10), |01) + |22), |02) 4- |20) and |11) + |20). The first one is the only 
product state in the kernel. Its multiplicity must be 6. 

For fc = 2 we give an example for each of the patterns (5, 1), (4, 2) and (3, 3). For the first example we set 

|V'o) = |12), |Vi) = |01)-|20), |V^2) = |02)-|21), |V'3) = |10)-|22). (47) 

The kernel is spanned by |00), |11), |01) + |20), |02) + |21) and |10) + |22). The first two of them are the only product 
states in the kernel. Their respective multiplicities are 5 and 1. 
For the second example we set 

|Vo) = |01)-|12), iV^i) = |10) - |21), |V^2) = |02)-|20), |7^3)-|22). (48) 

One can readily verify that kerp is spanned by |00), |11), |01) + |12), |10) + |21) and |02) + |20), and that |00) and 
1 11) arc the only product states in the kernel. Their multiplicities are 2 and 4, respectively. 
For the third example we set 

|Vo) = |02), |Vi) = |01)-|10), |V^2) = 111) - |20), |V'3) = |12)-|21). (49) 
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The kernel is spanned by |00), |22), |01) + |10), |11) + |20) and |12) + |21). The first two of them are the only product 
states in the kernel. Each of the two multiplicities is 3. 

For A: = 3 we give examples for each of the patterns (4, 1, 1), (3, 2, 1) and (2, 2, 2). For the first example we set 



|Vo)-|02), |V^i)-|20), |V^2)-|01)-|12), |7^3) = |10)-|21). 



(50) 



The kernel is spanned by |00), |11), |22), |01) + |12) and |10) + |21). The first three of them are the only product 
states in the kernel. Their respective multiplicities are 1,4,1. 
For the second example we set 



|V^o)-|20), iV^i) - |01)-|22), |V^2) = |02)-|12), j^^s) = |10) - |21). 



(51) 



The kernel is spanned by |00), |11), |02) + |12), |01) + |22) and |10) + |21). The first three of these pure states are the 
only product states in the kernel. Their respective multiplicities are 1,2 and 3. 
For the third example we set 



H'o) = \n), |V'i) = |02)-|12), |V'2) = |20)-|21), It^s) = |01) + |10) + |22). 



(52) 



The kernel is spanned by |00), |02) + |12), |20) + |21), |01) + |22) and |10) + |22). The first three of these pure states 
are the only product states in the kernel. For each of them the intersection multiplicity is 2. 

For A: = 4 we give examples with intersection patterns (3, 1, 1, 1) and (2, 2, 1, 1). For the first example we set 



|^o) = |01), |V'i) = |G2)-|ll) + |20), IV-2) = |10)-|22), l^^s) = |12) - |21). 



(53) 



The kernel is spanned by |00), |02) + |11), |11) + |20), |10) + |22) and |12) + |21). The first product state in the kernel 
is 1 00) and the other three are given by rank one matrices: 





1 C 

c e 1 



For 1 00) the intersection multiplicity is 3, and for the other three points it is 1. 
For the second example we set 



(54) 



1^0) -|10), iV^i) = |01)-|22), |V^2) = |02)-|12), |7^3) = |20)-|21). 



(55) 



The kernel is spanned by |00), |11), |02) + |12), |20) + |21) and |01) + |22) . The first four of these pure states are the 
only product states in the kernel. Their intersection multiplicities are 1,1,2 and 2, respectively. 
For k = 5 there is only one possible intersection pattern, namely (2, 1, 1, 1, 1). We set 



|^o) = |01)-|20), |Vi) = |02)-|ll), |V^2) = |10)-|22), |7^3)-|12)- 

The kernel is spanned by the states |00), |01) + |20), |02) + |11), |10) + |22) and |12) + |21) 
in the kernel is |00) and the other four are given by rank one matrices: 



|21). (56) 
The first product state 

(57) 



where i is the imaginary unit. For |00) the intersection multiplicity is 2, and for the remaining four product states 
each multiplicity is 1. 

For k = 6 there is again only one possible intersection pattern, namely (1,1, 1,1, 1,1). We set 



" 1 


1 


1 " 




1 


1 


-1 " 




i 


—i 


1 ' 




—i 


i 


1 ' 


1 


1 


1 




-1 


-1 


1 




-1 


1 


i 




-1 


1 


—i 


1 


1 


1 




1 


1 


-1 




—i 


i 


-1 




i 


—i 


-1 



|^o)-|01)-|12), iV^i) |02)-|21), |V^2) = |10)-|22), l^^s) - |11) ^ 

The kernel is spanned by the states |00), |01) + |20), |02) + |11), |10) + |22) and |12) + |21) 
in the kernel is |00) and the other five are given by rank one matrices: 

1 ^ 

e e e 
e c e 



|20). (58) 
. The first product state 

(59) 
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Clearly, in this case each multiplicity must be 1. □ 
The above example for k = b shows that a state of rank 4 whose range contains no product state may fail to be a 

PPTES. Indeed, the kernel of a PPTES of rank 4 contains exactly 6 product states (see Theorem [22] below) . 

On the other hand the example that we chose for fc = 6 is SLOCC-equivalent to the Pyramid UPB in [l^l- One way 

to see this is simply to verify that the quintuple of product states given by the above matrices for ^ = exp(27rzA:/5), 

/c = 0, . . . , 4 and the UPB quintuple for the Pyramid (see Eq. (|64l) below) have the same invariants. If we exchange 

the parties A, B and transform the product states in the kernel by the local operator 5* (gjS', 5 = \ ® \/l + VS 

we again obtain the Pyramid. 

Of course, the intersection of and Yj2,2 does not have to be proper, i.e., it may have positive dimension. For the 
comparison with Theorem 1101 we provide a scenario where there are infinitely many product states in the kernel. 

Lemma 11 Let p be a 3 x 3 state of rank 4 such that, rank(a;|p|a;) = 1 for some \x) G Ha- Then ker p contains 
infinitely many product states. 

Proof. Since the operator <t = {x\p\x) has rank one, it suffices to observe that the 2-dimcnsional subspace \x) 0kcr(<T) 
is contained in the kernel of p. □ 



B. Rank-4 PPTES with no product state in the range 



In this subsection we consider a related problem of describing the product states in the kernel of states of rank 4 
having no product state in the range. The set of such states properly contains all 3 x 3 PPTES of rank 4, as one will 
sec later in Theorem 1251 It is thus important to have a general understanding of this set. 

Lemma 12 Let p be a 3 x 3 state of rank 4- Then TZ{p) contains at least one product state when its kernel contains 
either 

(a) two linearly independent product states \a,b) and \c,d) with \a) = |c) or \b) = \d) , or 

(h) three linearly independent product states \ui,Vi), z = 1,2,3 such that the \ui) or the \vi) are linearly dependent. 

Proof. If (a) holds, say \a) = \c). then the 7-dimcnsional space span{|a, &), |a, d)}-'- contains Tl{p) and the 6- 
dimensional subspace V ^\a)^ ®T-Lb. Since V {^TZ{p) has dimension > 3, it must contain a product state. 

If (b) holds, say the \ui) are linearly dependent, then the 6-dimensional space {\ui,Vi) : i = 1,2,3}"'" contains TZ{p) 
and the 3-dimensional subspace V = {\ui),i = 1,2,3}""" ^Hb- Hence, the subspace V r\TZ{p) has dimension > 1. 
Since each nonzero vector in y is a product state the assertion follows. This completes the proof. □ 

Note that the converse does not hold; i.e., both (a) and (b) may fail even though 7?.(p) contains a product state 
(see the first example for k = 3). Nevertheless, this result will be strengthened in the case of PPTES in Sec. [V] 

On the other hand there may exist only 3 product states in general position in a 5-dimensional kernel, when there 
is no product state in the range of p. We consider 5 states in kerp represented by the matrices Ci,i = 0, 1, 2, 3, 4 



(60) 



respectively, where a 7^ 0,6 ^ 1. Up to ILOs, this is also the generic expressions of basis in a 5-dimensional kernel 
where there are only 3 product states in general position and there are no product states in the range of p. 

Lemma 13 For any 3x3 state p of rank 4, whose kernel is spanned by the pure states i6(J\) . TZ{p) contains no product 
state while ker p contains only 3 product states and they are in general position. 

Proof. First we show that there is no product state |a, b) e TZip). Because the first three product states in the kernel 
are |00), |11), |22), we have that \a,b) = \i){xi\{i + l) mod 3) + yi\{i + 2) mod 3)), or {xi\{i + I) mod 3)+yi\{i + 2) 
mod 3))\i), i = 0, 1,2. Since one of them is orthogonal to the latter two states 6*3,(74 in Eq. (pO)) . there are two 
parallel rows or columns in the same position of 6*3,(74. This is evidently impossible, so there is no product state in 

n{p). 

Second wc show that |00), |11), |22) are the only three product states in the kernel. It is easy to see that wc must 
use C4. Wc compute the linear combination oi Ci,i = 0, 1, 2, 3, 4 such that 



" 1 





0" 




" 





0" 




' 





0" 




" 





1 " 




" 


a 


b 
















1 
























1 




a 





-1 + 6 






























1 




1 


1 







1 









C := uCo + vCi + WC2 + xCs + C4 



u a b + X 
a V — 1 + 6 + a; 
1 + X X w 



(61) 
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which has rank 1 when it is a product state. Evidently x ^ —1,0. Then we can deduce that v = which leads to 



det 



a b + X 

V —1 + b + x 

Furthermore we consider 5 states in ker p represented by the matrices C^, i = 0, 1, 2, 3, 4 



7^ 0. Hence C cannot be a product state Vw, v, w, x. This completes the proof. 



□ 



" 1 





o" 




" 





o" 




" 


1 


1 " 




" 


1 


1 




" 


'-^ 2+c 


c 
















1 













1 















1 




































1 




1 


-2-c 


-1-c 







4c 
2+c 


c 



(62) 



respectively, where the real c ^ —2,0. Then we have 



Lemma 14 For any 3x3 state p of rank 4, whose kernel is spanned by the pure states \6'^) . TZ{p) contains no product 
states while kerp contains only 2 product states (up to scalar multiple) and these two states are in general position. 

Proof. The proof is similar to that for Lemma [T51 First we show that there is no product state \a,b) € Ti-ip)- 
Because the first two product states in the kernel are |00), |11) we have that |a, b) ~ (a;o|0) + X2|2))(?/i|l) + y2|2)), or 
{xi\l) + a;2|2))(yo|0) + J/2|2)), or |2, b) or |a, 2). One can check that none of them exists, so there is no product state 
in Tl{p). Second we show that |00), |11) are the only 2 product states in the kernel. To simplify the proof, we notice 
that the the first 4 blocks cannot generate new product states. So we must need C4. Further we consider three cases 
for the linear combination C := uCq + vCi + wC-z + xCs + C4, namely g = 4c/(2 + c)^, g = 1 and the rest. In each of 
these cases, one can easily show that C cannot be a product state. This completes the proof. □ 
We lack examples of states whose range contains no product state, and whose kernel contains only 1 product state. 
This problem, as well as Lemmas [T3] and [T31 is more relevant for the characterization of NPT states, which is an 
essentially useful quantum- information resource (for a recent paper see 0). From the next section we will focus on 
the main topic of this paper, namely the description of the 3x3 PPTES of rank 4 via the UPB construction 0]. 



IV. CHARACTERIZATION OF EQUIVALENT 3x3 UPB 

Let us denote by £4 the set of all PPTES of rank 4 in a 3 (g) 3 system. Our main objective is to prove a conjecture 
which was raised in [3l| and gives a full description of the set £4. This has close connection with the family of PPTES 
constructed in [l3| via UPBs. 



A. PPTES of rank 4 and UPB 



Let hi denote the set of all UPBs in H = Ha ^ Hb- We denote by W the set of all quintuples (tp) := (IV'fc) = 
\4'k) <8i \Xk))k=o of (normalized) product states such that the set {|V'fc) : < fc < 4} is a UPB and the following 
orthogonality relations hold: 

(0»l0j+i) = {Xi\Xi+2) = 0, (63) 

where the indexes are taken modulo 5. We have a natural projection map U° — > U which associates to a quintuple 
(-0) = (|-0fc))t=o ^ ^° the UPB {V'} := {\ipk) : < k < 4} eU. It was shown in ^ that this map is onto. It is 
clearly 10-to-l map because the cyclic permutation of the {ipi) and also the reflection which interchanges the indexes 
via the permutation (0)(14)(23) has no effect on the set {ip}, and leaves U° globally invariant. 

There is a natural map II : U —?' £4 which associates to G U the state n{?/'} := (1/4)P, where P is the 
orthogonal projector of rank 4 with kerP = span{'0}. The following conjecture, which gives explicit description of £4 
was raised recently in [3T| and supported by vast numerical evidence. 

Conjecture 15 Every state p £ £4 is the normalization of A®B n{-0} A"^ ® B'^ for some {A,B) E GL and e U. 
The proof of this conjecture will be given in Theorem [25] of Section [V] 
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We fix o.n. bases {|0)a, |1)a, |2)a} and {|0)_b, |1)b, |2)b} of T-La and "Hs, respectively. It was shown in [l^ that 
every quintuple \'4'k) = l^fc) ® \Pk)^ fc = 0, . . . , 4, in 11° is LU-equivalcnt to one in the following 6-paranieter family: 

|ao> = |0>A, (64) 

|ai> - \1)a, (65) 

\a2) = cosOa^a + sin0^|2)^, (66) 

lag) = sin7Asin6'^|0)A - sin7ACOs6'A|2)A +cos7Ae"^'^|l)^, (67) 

|a4) = (sin7ACOS0Ae"^-*|l)A + cos7A|2)A) , (68) 
i\a 

l/3o) - (69) 

|/3i) =sin7Bsin0B|O)B-sin7BCOs0B|2)B+cos7Be^^«|l)B, (70) 

1/32) = |0)b, (71) 

I/33) =cos0B|O)B + sin0B|2)B, (72) 

I/34) = ^ (sin7BCOs0Be'^^|l)B +cos7b|2)b) ■ (73) 

iVB 

The 6 real parameters are the angles: 7a, s, Oa,b and c/'a,b, and the normalization constants Na,b are given by the 
formulae 



Nam = Y cos2 7A,B + sin 7a,b cos^ Oa,b- (74) 

The first four angles are required to have nonzero sine and cosine, while the angles 4>a,b niay be arbitrary. It is not 
hard to show that the parameter domain can be further restricted as in the following lemma. 

Lemma 16 Every quintuple (ip) = (|^fc))|.=o ^ iV'fc) = \ctk) <Xi |/3fc) is LU-equivalent to one belonging to the 

family \64^ such that the four angles Ja,Bj Sa,b belong to the interval (0, 7r/2). 

Proof. We may assume that |ao) = |0)a, \cii) — |/3o) = |l)s and |/32) = |0)b- As \a2) is a linear combination 

of \Q) A and |2)a, we can choose the overall phase of \a2) so that the coefficient of \Q) a is positive. By applying 
a diagonal unitary matrix Ua = diag(l,l,^), we can also assume that the coefficient of \2)a is positive. Thus 
\a2) = cos 6*^10) A + sin 6*^12) A for some 9a € (0,7r/2). 

Next, jaa) is a linear combination of |1)a and sin6'A|0)A — cos6'a|2)a- Wc can choose the overall phase of {a^} so 
that the coefficient of sin 0a |0) a — cos6'a|2)a is positive. Thus 

lag) =sin7Asin6'^|0)A-sin7ACOs6lA|2)A + cos7Ae*'^'*|l)A (75) 

for some 7a € (0,7r/2) and some angle 0a- 

Finally, [04) is a linear combination of |1)a and |2)a. Wc can choose the overall phase of \a2) so that the coefficient 
of |2)a is positive. Since |q!4) is orthogonal to jaa), we have 

l"4) = (sin7ACOS0Ae"^-*|l)A + cos7A|2)A) (76) 

iVA 

with 7a G (0, 7r/2) and the positive normalization constant Na- 

The same arguments can be used on Bob's side. □ 
We shall give just one example, namely the UPB quintuple known as Tiles. Its parameters, as given in [13, p. 

394], are (t>A.B = 0, Oa.b = 7a, b = 37r/4. This quintuple is LU-equivalent to the one given by parameters (f>A,B = 

9a, B = 1A,B = 7r/4. The local unitary transformation that we can use in this case fixes the vectors |0)a, |1)a and 

|0)b, |1)b, and sends |2)a and |2)b to their negatives. 
Let F be the subfamily of the family ([M)) obtained by restricting the domain of parameters so that the four angles 

7a, b, 9a,b belong to the interval (0,7r/2) while the angles 4'a,b belong to (— tt, tt]. Since the domain of parameters 

is connected, the family F is also connected. By Lemma [12] we have U° = U(3) x U(3) • J-, and so the set U° is 

connected too. 

Lemma 17 Assume that A(E)B{ip) = where {A, B) € GL and the quintuples (ip) = {\ipk))t=o '^"■'^ (^') ~ (IV'fc))fc=o 
belong to IA° . Then there is a positive constant c such that cA and c~^B are unitary. In particular, if (-0) and (ip') 
are SLOCC-equivalent then they are LU-equivalent. 
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Proof. Let us write 

Hk) = \<i^K)®\Xk), Wk) = Wk)®Wk). fc = 0,...,4. (77) 

By the hypothesis we have 

A\^k)<^B\xk) = Wu)^\x'k). fc = 0,...,4. (78) 

Hence, A\4>u) = au\4>'f,) and B\xk) ~ o.k^\x'k) some scalars a^. Since {4>'i:\4>'k+i) ~ 0' deduce that 
(0fe|A'''A|(/)fc_|_i) ~ 0. Consequently, A must map the plane spanned by |(7!>fc+i) and |0a;-i) onto itself. This clearly 
implies that A is a scalar matrix, i.e., there is a scalar c > such that cA is a unitary matrix. A similar argument 
shows that c~^B is also unitary. □ 
Note that SLOCC-equivalence is different from the BP-equivalence which does not require the identical global scalar 
for simultaneous transformations A (g) Blipk) = iV'fe): fc = 0, • • • ,4. By Lemma ITTl two SLOCC-equivalent quintuples 
(-0) and (V-'') in T must be connected by local unitary operators Ua <E} Ub- 

Proposition 18 // two quintuples {\tpk)), ilfpk)) £ ^ '^''e SLOCC-equivalent, then \^k) — iV'fe) for fc = 0, . . . , 4. 

Proof. As in Eqs. (|M)) wc write \ipk} = \ak) ® |/?fe) and similarly let \tp',.) ~ jaj.) ® |/3^). Let {ja.b,9a,b,4'a.b) and 
ilA Bt^'a 'I^'a b) be the parameters of the two quintuples, respectively. 

The SLOCC-equivalence implies that the quintuples {\oik))k={) and (|q!/;))|=o projectively equivalent. We point 
out that these quintuples are in general position and so, by Proposition [3] and O they must have the same invariants 
J^. By using the expressions in Eq. ([64]) and the formulae Eq. ([2]), we find that these invariants for [\'4>k)) are given 

by 

Ji ^- tan' cos' 9 a, J2 = coie^ ' -^3^ = ^ cot' 7^, 

J? = -4«-, J? = Ta\" 2 , ^1 + cos' 9b tan^ 75 

-L sin'^ ^ 1+cos'^ tan^ 7B ' ^ -d /r> 

Similar formulae are valid for the quintuple (1^^))- The equalities J/^ = and — imply that 

cos' U = cos' , C = 7, 0; X = A, B. (80) 

Since all these angles belong to (0, 7r/2), wc conclude that ^'a b — 1a,b and 9'^^ g — (^a,b- It remains to prove that 
also (f)'j^ B = (t^A^B- 

By Lemma ITTl there exist unitary matrices A = [aij] and B = [bij] such that 

A\ak) ® i?|/3fe> - 14) ® 1/3^) =0, fc = 0, . . . , 4. 

We view each of these equations as a matrix equation. For fc = it gives 011622 = 1 and 021 = 03 
Thus we may assume that an = 622 = 1 and, since A and B arc unitary, we must also have ai2 = ai 
For fc = 2, the matrix equation shows that 611 = 1 and 033 = 1. We deduce that 023 — 032 — 613 
Now, for fc = 1 the matrix equation implies that 022 ~ 633 = 1. This means that A = B = It, anc 





(81) 


I — bi2 — 632 


= 0. 


3 = 621 = ^23 


= 0. 


= 631 = 0. 




1 so 4)'yi B 4 


'A,B- 




a 



B. UPB symbols 

Let us say that a 5-dimensional subspace C H is of UPB type if it is BP-equivalent to a subspace spanned by a 
UPB. We can characterize the UPB-type subspaces by using invariants. For this purpose we attach a 6-letter symbol 
made up of letters N, P and p to each quintuple of product states in general position having all invariants real. 

Let us denote the open intervals {—00, 0), (0, 1), (1, +00) by the letters A^, p, P, respectively. {N is for "negative", 
p for "positive and small" and P for "positive and large".) Let (ip) := {\ipk))k=oi be any quintuple of product states 
in general position. For convenience, we shall say that (ip) is real if all of its invariants are real numbers. Since 
the invariants do not take values and 1, if (ip) is real its invariants must take the values in one of the intervals 
N,p,P. For real (■0) we define its symbol to be the 6-letter sequence obtained by replacing each of its invariants 
J^, J,^, J^, Ji , J2 , J3 by the interval to which it belongs. For instance, each quintuple (tp) <S has NPNPpP as 
its symbol because <0, > 1, J3 < 0, Jf > 1, < < 1 and > 1. Altogether there are 144 symbols that 
arise in this manner. Each symbol can be broken into two parts, A and B. The A [B] part consists of the first [last] 
three letters. Because of the identity J1J2J3 = 1, there are only 12 possibilites for each of the two parts: 

NNP, NNp, NPN, NpN, pNN, pPP, pPp, ppP, PNN, PPp, PpP, Ppp. (82) 

We shall refer to the 12 symbols in Table 1 as the UPB symbols. 
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Tabic 1: UPB symbols and associated permutations 



NNPPPp (12) (34) 
PNNpPP (01) (34) 
pNNPpp (01) 



NNpppP (12) 
PPpNNP (23) 
pPPPNN (13) 



NPNPpP id 
PpPNPN (1243) 
pPpNpN (03) (24) 



NpNpPp (34) 
PpppNN (02) (14) 
ppPNNp (01) (23) 



We can now prove the main result of this section. 

Theorem 19 Let W QH be a 5- dimensional subspace containing exactly 6 product states (up to scalar factors) and 
assume that these states are in general position. Denote by ^> the collection of the 720 quintuples (ijj) := (|'0fc))fc=o' 
selected from these 6 product states. 

(a) Assume (ijj) £ ^/ is real and its symbol is UPB. If a is the permutation from Table 1 associated to this symbol, 
then i\i^ak))t=o symbol NPNPpP. 

(b) If some (ip) € is real and its symbol is UPB, then W has UPB type. 

(c) Conversely, ifW has UPB type, then all (ip) € ^ are real and their symbols are UPB. 

Proof. The assertion (a) can be proved by straightforward verification. We shall give details for one case which we 
shall need later. Assume that {tp) has symbol pNNPpp. Since we work with non-normalized states, without any loss 
of generality we may assume that (ip) is given by the pair of matrices 



A 



10 11 

1 1 -a 
11-6 



B 



10 1 
10 1 
11 



1 

c 
1/d 



By computing the invariants of (ip) we obtain that 



= a/b, = -b, = -l/a, Jf = cd, Ji = l/d, jf = 1/c, 



(83) 



(84) 



and so we must have 6 > a > and c,d> 1. From Table 1 we find that a = (01) in this case. Hence, the quintuple 
(-0') := {\'ipak))k=o given by the matrices A' and B' which are obtained from A and B, respectively, by switching 
the first two columns. For the invariants of (ip') we obtain the formulae 



= -1/6, J^-" = b/a, 4^ = -a, = d, = 1/cd, 



(85) 



Thus, the symbol of (V/) is indeed NPNPpP. 

Next we prove (b). By using (a) we may assume that {ip) has symbol NPNPpP. Since JfJ^J^ = 1 and 
Ji J2 Jz = I7 the equations ([7^ can be solved for the four angles 'yA,dA,'^B,(^B- Now the assertion follows from 
Propositions [3] and [5j 



To prove (c) we observe that any pair of real matrices 





U 



1 a 13 
10 1a 
^ -a 1 



V 



1 (5 7 
1 1 7 
-7 1 J 



(86) 



with nonzero parameters a, /?, 7, S defines a non-normalized UPB given by the tensor products of the corresponding 
columns of U and V. Moreover, any UPB is BP-equivalent to one of this form. The invariants of the above UPB are 



1 



1 



5^ + 



(87) 



Evidently, this quintuple has the symbol NPNPpP. By using these invariants and the formulae Eqs. ^ and 
from Proposition [7] we find the sixth product state in W and extend the matrices U and F to 3 x 6 matrices 



U = 



V 



1 a 13 a[{l + a^){l + j^ +5'^)+ f3^j^ +6^)]/[l3{l + a^ 

10 1 a a(l + 72 + 52)[^2 ^ (^,2 ^ ^2)(^2 ^ j2)]/[^2^2 ^ (^2 ^ J2)(^2 ^ ^2(^2 ^ ^2))J 

^ -a 1 1 

15 7 7[a'(l + a' + /3')(l+7' + '52)+/3'(l + 7')]/[/3''5(l + a2+72)] 

1 1 7 ^{a'+l3^)[{l + a^){l+j' + S^)+l3^{j^ + 5^)]/[/3^S^ + {a^+l3^)i-r^+6^))] 

-7 1 (5 1 



(89) 



16 



by appending this new product state. 

The symmetric group 5*6 permutes the 6 product states and induces a permutation representation on the 720 
quintuples made up from these 6 product states. For instance, if we choose the quintuple corresponding to column 
numbers 3,6,2,1,5 (in that order) then the invariants are 

jA ^ (l + «^+7^)(l + 7^ + <^^) 

^ (l+72)[(l+a2)(l+^2 + 52)+^2(^2+^2)]' 

J2 = T^i^, (91) 



1 + 7^ + ^ 

(1 + a2)(l +72+(52)+ /?2(72 + ^2) 

1 + a2 + 

/32(72+j2)(i^^2^^2) 



"'s - 1 , „,2 1 ^,2 ' vJ^; 



a2^2[(l + ^2)(l+^2 + J2)+^2(^2+^2)]' (93) 

^2^ ^ -f^!(i±^, (94) 

,B ^ 7^[(l + C.2)(l+72+^2)^^2(^2^^2)] 

(7' + '52)(l + a' +/3')(l + a2 + 72) ^ ^ 

and the associated symbol is ppPNNp. 

A brute force computation shows that there are only 12 different symbols that belong to these 720 quintuples, 
namely the UPB symbols listed in Table 1. This completes the proof. □ 

It is easy to see that all 144 symbols arise from some real quintuples of product states. Hence, apart from PPTES, 
the NPT states may also have 5 dimensional kernels with exactly 6 product states in general position. 

By means of Theorem [HI we can operationally decide the UPB-type of 5-dimensional subspaces. This is the key 
tool in our proof of Conjecture [15] in the next section. 

V. DESCRIPTION OF 3 x 3 PPT STATES OF RANK 4 

We shall first analyze the 5-dimensional subspaces which arise as kernels of 3 x 3 PPT states of rank 4. This will 
help us to resolve a conjecture proposed in [sij . 

A. Product states in the kernel of 3 x 3 PPT states of rank 4 

We need the following lemma which we proved recently in 0, Lemma 20]. 

Lemma 20 Let p be a 3 x N state such that for some \a) € Ha, rank(a|p|a) = 1. If p is NPT then it is distillable. 
If p is PPT and N = 3, then p is separable. 

Let us first handle the separable states of rank 4. 

Lemma 21 Let p = X]^=o 1^*' bi){ai, bi\ be a separable 3x3 state of rank 4- If the four product states \ai, hi) are not 
in general position, then kcrp contains a 2- dimensional subspace V ®W with V C "Ha and W C'Hb (which consists 
of product states). Otherwise (a) kerp contains exactly 6 product states, and (b) these 6 product states are not in 
general position. 

Proof. Assume that the \ai,hi) are not in general position. We consider first the case where two of the \ai) or two 
of the \bi) are parallel, say |ao) cx |ai). If \y) £ Hb is orthogonal to I62) and I63), then |ao)''" (8) \y) C kerp. In the 
remaining case we may assume that, say |ao), |ai) and [02) are linarly independent while I62) belongs to the span of 
|5o) and If the state \y) S Hb is orthogonal to |6o) and |&i), then las)-'- (X) \y) C kerp and so the first assertion is 
proved. 

Next assume that the \ai,bi) are in general position. By the Four Point Lemma, we may assume that these product 
states are in the canonical form, i.e., \ai,bi) oc \ii), for i = 0,1,2 and [04,64) = X]?j=o I*-?)- "^^"^ ^^-^ product states 
K) ^ 1^')) and (|j) — |fc)) (g) |z), with («, j, k) a cyclic permutation of (0, 1, 2), belong to kerp. To prove (a) one 

just needs to verify that there are no additional product states in ker p. We omit the details of this verification. For 
(b), it suffices to observe taht the sum of |0) — |1), |1) — |2) and |2) — |0) is 0. This completes the proof. □ 
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Prom now on we focus on the PPTES of rank 4. We recall from 01 that any state p of rank 4 acting on a 3 (g) 3 
system can be written as 

2 

P= E I*)(J'I®C'1C,' (96) 

where the blocks Q are 4x3 matrices. 

Theorem 22 The kernel of any 3x3 PPTES of rank 4 contains exactly six product states. Moreover, these six states 
are in general position. 

Proof. Let p be a PPTES of rank 4. Assume that ker p contains infinitely many product states. By Lemma [12] they 
are all in general position. Hence we can apply Proposition [7] to any quintuple of product states in ker p. As ker p 
contains infinitely many product states, only the case (a) of that proposition applies. The third assertion of that case 
contradicts our hypothesis that p is PPT. This contradiction shows that kerp contains only finitely many product 
states. 

We may assume that p is written as in Eq. (|96p . We have p = C'^C, where C ~ [Cq Ci C2]. We shall simplify this 
expression by using the techniques similar to those in . We can replace C with U C where [/ is a unitary matrix, 
without changing p. The effect of a local transformation p ^ {I ® B)'^ p {I ® B) \s to replace each Ci by GiB. 
Similarly, a local transformation p — > (A (g) p {A® I) acts on p via block-wise linear operations, such as adding a 
linear combination of Co and Ci to C2, etc. We can apply these kind of transformations repeatedly as many times as 
needed. Note also that if the jth column of Ci is then \ij) £ ker p. By Lemma each block Ci has rank at least 2. 

Since kerp has dimension 5, it must contain a product state. We choose an arbitrary product state in ker p. 
By changing the o.n. bases of Ha and "Hs, we may assume that the chosen product state in kerp is |00). Since 
= p|00) = J2i N) ® C'|Co|0), we must have CjCo\0) = for each i. In particular, CqCoIO) = which implies that 
Co|0) = 0, i.e., the first column of Co must be 0. Hence, the block Co must have rank 2, and we may assume that 



Ci = C2 = [c„]. (97) 



Let (T = p^ and observe that its first entry is 0. Since a > 0, the first row of a must be 0. We deduce that 
bii ~ b-ii = cii = C21 = 0. Since ker p contains only finitely many product states, the first columns of Ci and C2 must 
be linearly independent. By using an ILO on system A, we may assume that 631 = C41 = 1 and 641 = C31 = 0. Thus 
we have 

&12 bi3 C12 

^ ^ &22 &23 C = ^ ^'^^ 

1 ^32 ^33 ' C32 
_ &41 643 J [1 C42 

while Co did not change. 

Let be the 4-dimensional complex projective space associated to ker p. The Segre variety E2.2 and this 
intersect properly and we claim that the intersection multiplicity at the point P = |00) is 1. To prove this claim, we 
introduce the homogeneous coordinates for the projective space associated to H: If \ip) = X]?j=o '^ijl'^j) then 
the homogeneous coordinates of the corresponding point G are ^ij ~ Oij. The computation will be carried 
out in the affine chart defined by ^00 7^ which contains the point P = |00). We introduce the affine coordinates 
Xij-i {h j) 7^ (O7O), in this affine chart by setting Xij = ^ij/^oo- Thus P is the origin, i.e., all of its affine coordinates 

Xij — 0. 

The computation of the intersection multiplicity is carried out in the local ring, say R, at the point P. This local 
ring consists of all rational functions f /g such that g docs not vanish at the origin, i.e., / and g are polynomials (with 
complex coefficients) in the 8 affine coordinates Xij and g has nonzero constant term. By expanding these rational 
functions in the Taylor series at the origin, one can view i? as a subring of the power series ring C[[a::ij]] in the 8 affine 
coordinates Xij . We denote by m the maximal ideal of R generated by all Xij . 

The range of p is the 4-dimensional subspace spanned by the pure states lipi), « = 1, . . . , 4, given by the four columns 
of C^ (see the proof of Proposition 6 in Q). In the matrix notation, these pure states are represented by the following 



Co 



1 

1 







Cl3 

C23 
C33 
C43 



(98) 
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matrices 



'o 


1 







'o 





1 




"o 






















bh 


h* 







h* 
"22 


h* 
"23 




1 


h* 
"32 


h* 
"33 







h* 
"42 


h* 
"43 





r* 


r* 







r* 
'^22 


r* 







r* 
•-32 


r* 

'-33 . 




1 


r* 
'-42 


r* 

•-43 . 



Since kerp = TZ{p)^ , the subspace V'^ is the zero set of the ideal Ii generated by the four Unear polynomials: 

xoi + &i2a;ii + bi3Xi2 + £122:21 + C13X22, 
X02 + &22a;ii + 6232^12 + C22a;2i + C23a;22, 

XlO + b^2Xll + 6332^12 + C32a;21 + C33a;22, 

X2n + 6422:11 + 6432^12 + C422;2i + 0432:22- 
The piece of the Segre variety contained in our affine chart consists of all matrices 

1 2:01 2:02 
xio 2:11 X12 
X20 X21 X22 

of rank 1. It is the zero set of the ideal I2 generated by the four polynomials: 

xii ~ xoixio, 2:12 - 2:02X10, 2:21 - 2:01X20, 2:22 - 2:022:20 

^2 



(99) 



(100) 
(101) 
(102) 
(103) 



(104) 



(105) 



The quotient space m/m is an 8-dimensional vector space with the images of the Xij as its basis. It is now easy to 
see that the images of the generators of Ii and I2 also span the space m/tn^. Hence, by Nakayama's Lemma (see 
p. 225]) we have Ii + I2 = na. Consequently, i?/(/i + 12) = C and so our claim is proved. 

Recall that we chose in the beginning an arbitrary product state in ker p and by changing the coordinates we were 
able to assume that this product state is |00). Since the intersection multiplicity is invariant under these coordinate 
changes, this means that we have shown that the intersection multiplicity is 1 at each intersection point of and 
^2,2- By Bezout Theorem the sum of the multiplicities at all intersection points is 6, and since all of the multiplicities 
are equal to 1 we conclude that the intersection consits of exactly 6 points. 

By Lemma [121 these six product states are in general position. This concludes the proof. □ 



B. r-invariant PPTES of rank 4 



We shall prove that every SLOCC-equivalence class in £4 contains a state which is F-invariant. 



Theorem 23 Any 3x3 PPTES p of rank 4 is SLOCC-equivalent to one which is invariant under partial transpose, 
i.e., for some {A, B) e GL and a = A®BpA'^®B'^ we have = a. 

Proof. By Theorem I22| we may assume that \ii) G kerp for i = 0, 1, 2. Hence we may assume that, in the formula 
(|96| for p, the column i + 1 of the block d vanishes for i = 0,1,2. By multiplying C = [Co Ci C2] by a unitary 
matrix on the left hand side and by performing an ILO with diagonal matrices we may assume that 



Co = 



" 


1 


h' 




"611 





bi3 




Cll 


C12 


0' 








1 


, Ci = 


621 





623 


, C2 = 


C21 


C22 
















631 





633 


C31 


C32 





_ 





0_ 




. 641 





bi3 _ 




. C41 


C42 


0_ 



(106) 



Let <T = p^ and observe that its first, fifth and ninth diagonal entries are 0. Since a > 0, the rows of a containing 
these entries must be 0. We deduce that 611, 621, cn, C21, 613 are and that the second column of C2 is orthogonal 
to the first and third columns of Ci. Since kerp contains only finitely many product states (up to scalar multiple), 
the first columns of Ci and C2 must be linearly independent. Hence, by applying a unitary transformation to the last 
two rows of the Ci and rescaling Ci and C2, we may assume that 631 = 0, 641 = 1 and C31 = 1. By the orthogonality 
property mentioned above, we conclude that C42 = 0. Thus we have 



Co 



' 


1 


b ' 




' 













C12 


0" 








1 


, Ci = 








&23 


, C2- 





C22 






















&33 




1 


C32 





_ 





0_ 




1 





643 . 




. C41 





0_ 



(107) 
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Since p is entagled and PPT, its range contains no product states. Consequently, the entries C12, C22, ^33, ^43 and 
h must be nonzero. We choose a phase factor zi such that €412:1 > 0. By niuhiplying the first columns of all Ci by 
zi and mutiplying Ci by zjf, we may assume that d :~ C41 > 0. By multiplying the last two columns of all Ci with 
1/^33 and multiplying Cq by 633, we may assume that 633 = 1. We choose the phase factor Z2 such that bz2 > 0. By 
multiplying the second columns of all Ci by Z2 and mutiplying Co by Z2, we may assume that 6 > 0. Since the last 
row of (T must vanish, we obtain that 643 = —l/d, C12 = — C22/& and C32 = — fe23C22- If ^23 = then also C32 = and 
\1lj4) — dlips) is a product state (see Eq. (|M|) '). Hence, c := 623 7^ 0. Next choose a phase factor Z2 such that cz2 > 0. 
By multiplying the last two columns of all Ci by Z2 and multiplying Ci by Z21 we may assume that c > 0. Finally, by 
multiplying the second columns of all Ci by a := I/C22 we have 



" 


a 


b' 




' 










' 


-1/b 0" 








1 










c 


, C2 = 





1 



































1 




1 


-c 








0_ 




1 





-l/d_ 




_d 


0_ 



(108) 



We claim that a has to be real. For that purpose we compute the principal minor of a obtained by deleting the first, 
fifth and ninth rows and columns. We obtain the expression c^(a — a*)^. Since this minor must be nonnegative, our 
claim is proved. 

It is now easy to verify that = p, which completes the proof. □ 
From the proof of Theorem E51 it follows that any PPTES p is SLOCC-equivalent to one given by Eq. ([M|) where 
the blocks Ci are given by Eq. (jlOSp . Moreover the parameters a,b,c,d are nonzero real numbers with b,c and d 
positive. The converse also holds (and it is easy to verify), i.e., if p is given by this recipe then it must be a PPTES. 
The above formulae (|108|) for the blocks C,; play an essential role in the proof of the next two theorems. 



Theorem 24 // p, p' € £4 have the same range, then p ^ p' . 

Proof. Clearly, in order to prove the theorem we can simultaneously transform p and p' by the same ILO. Thus we 
can assume that p is given by Eq. ([TO]) and that the blocks Ci are as in Eq. (|108p . For convenience, we shall not 
normalize neither p nor p' , and so we have to prove that p' is a scalar multiple of p. The range of p is spanned by the 
four pure states represented by the four matrices 






a b 










1 



































c 







1 




1 





-l/d 





-1/b 







1 







1 


-c 




d 









respectively. The first of these matrices is made up in the obvious manner from the first rows of the blocks Ci , and 
the other three matrices are constructed similarly. All this information is encapsulated in the matrix C = [Co Ci C2]. 

We also have p' = X]^=o where \ip'i) are four linearly independent linear combinations of the pure states 

Hence, the matrix C corresponding to p' is given by C = SC, where S is an invertible matrix of order 4. Thus we 
have p' = C'^HC where H = S'^ S is a positive definite matrix of order 4. If cr' = {p')^ , an easy computation shows 
that the first, fifth and ninth diagonal entries of a' are 0. Consequently all entries of a' in the first and fifth row must 
be equal to 0. These equations give immediately that H is a. scalar matrix, which completes the proof. □ 



C. Main result 



We can now prove our main result, i.e. Conjecture 1151 

Theorem 25 Up to normalization, any state p £ £4 has the form A(E) B (E) for some [A, B) G GL and 

some {■!/;} e Li . 

Proof. The proof is based on Theorem [THl The first step is to construct the 6 product states in ker p. Like in 
Theorem [24l we assume that p is given by Eq. (|96)) and that the blocks Ci are as in Eq. (|108p . A direct computation 
shows that the 5-dimensional space kerp is spanned by |00), |11), |22), and any two of the pure product states 



((1 + 52 + 62^2)^^ _ ^2^) 1^^^^ ((^ ^ ^2 ^ ^2^2)^^ _ ^2^-) C - (1 + C^)z, 

(czi — 1 — d'^) / d Zi{czi — 1 — d^)ld czi — 1 

1 Zi d{czi — \) / {czi ~ 1 — d^) 



(110) 
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where Zi (i = 1, 2, 3), are the three roots of 

abz{cz - 1 - d'^){c - (1 + c^)z) + d{cz - l){h'^c -{l + b^ + l?c^)z) = 0. (Ill) 
The l.h.s. is a cubic polynomial, say /(z), in the unknown z. We have 

/(O) = -h^cd < 0, (112) 

1 + c/ [1 + c^)^ 

/(i±^) = . (62^2 _ (1 + 52 + b^c^){l + rf2)) < 0. (114) 
c c 

Since < c/(l + c^) < (1 + d^)/c, we deduce that Eq. (|llip has three distinct nonzero real roots. One of them is 
zi e (0,c/(l + c^)), the second Z2 £ (c/(l + c^), (1 + d'^)/c), and the third is Z3 < if a > and Z3 > (1 + d'^)/c if 
a < 0. We also have 



/(1/c) = a6dVc2 ^ 0, (115) 
-b^){l + d^)+b^ 
(1 + 62 + &2c2)3 



,3 2 (1 + 6^)(l + rf^) + 6W ^ 



where A = 6^c/(l + 6^ + 6^c^). Hence, for a > we have 

Z3<0, A < zi < c/(l + c2), l/c< Z2 < (l + rf^)/c, (117) 



and for a < 



0<Zi<A, c/(l + c2) < Z2 < 1/C< (l+d2)/c< Z3. (118) 



The second step is to compute the invariants for one of the quintuples selected from the above 6 product states. 
For that purpose we shall use the quintuple (|00), |11), |22), \ip2))- A computation gives the formulae: 

tA _ 1 + - CZ2 _ Z2{zi - A) „ ^1(1 + - CZl)(z2 " A) 

" l + d2_czi' " Z1(Z2-A)' - ^^(^^^2_^^^)(^^_^)' i^^^i 

.5 ^ 2:2(1 - CZ2)(CZ1 - 1) B ^ {I + d^ ~ CZI){CZ2 ~ I) jB ^ ^ /ir,r|N 

"'I Zi(l+d2_cZi)(cZ2-l)' ' (l+d2_cZ2)(cZi-l)' ^ ' ^^^^^ 

As the third step, we have to compute the symbol associated to this quintuple. There are two cases to consider 
according to whether a > or a < 0. We claim that the symbol is ppPNNp in the former case and pNNPpp in the 
latter case. This verification is of routine nature and Figure 1 may be useful for that purpose. We shall just verify 
the claim in the case a > 0. In that case, it is obvious that < < 1, < 0, < 0, and < < 1. Since 
A < zi < Z2 and the function t/{t — A) is decreasing for t > A, we indeed have < < 1. Since JiJ2Jz ~ we 
deduce that > 1. Thus we have shown that the symbol is ppPNNp if a > 0. 
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Finally, since ppPNNp and pNNPpp are UPB symbols, wc conclude that ker p is a S-dimensional subspace of UPB 
type. Hence, we can now apply Theorem [Ml to complete the proof. □ 

It follows from the above proof that every 3x3 PPTES of rank 4 is SLOCC-equivalcnt to one given by Eqs. (|96|) 
and (jlOSp with positive a, 6, c, d. 

The next corollary shows that there is no way to single out one of the six product states in the kernel of a PPTES 
of rank 4 in the sense that any quintuple of these states is BP-equivalent to a qunituple formed from the five states 
of a UPB. 



Corollary 26 Let p G £4 and let \tpk), k ~ 0, . . . ,A, 
(A, B) G GL such that the product states A® B \i}}k) 



he any five of the six product states in ker p. 
fc = 0, . . . , 4, form a non-normalized UPB. 



Then there exists 



Proof. By Theorem[251 ker p is a 5-dimcnsional subspace of UPB type. By Theorem[Tni the symbol of the quintuple 
dV')) = (IV'fc))fc=o ^ UPB symbol and there exists a permutation (fco, . . . , ki) of the indexes (0, . . . , 4) such that the 
symbol of the quintuple (IV"')) (l'0fci))f=o NPNPpP. We can now conclude the proof by using the argument 
from the last paragraph of the proof of Theorem [T9l □ 
We now analyze the stabilizer oi p £ £4 in the product PGL = PGL3 x PGL3 of two projective general linear 
groups. Thus we have to consider {A, B) E GL such that A ® B p A^ ® B'^ ^ cp ioi some scalar c > 0. 

Proposition 27 The stabilizer Gp of any p £ £4 in PGL is a finite group isomorphic to a subgroup of the symmetric 
group Sq. In the generic case the stabilizer is trivial. 

Proof. Let us denote by P^, i = 1, . . . , 6, the six points in the projective space associated with ker p. Assume 
that {A,B) £ GL maps p to cp for some c > 0. Then A® B must leave invariant Ti{p) and {A~^ (E) B~^)'^ must 
leave invariant kerp and permute the 6 product states, i.e., the points Pi. The map which assigns to {A,B) this 
permutation, say tta.b € 5*6, is a group homomorphism. If tta.b is the identity permutation, then the Four Point 
Lemma implies that {A~^ (g) B~^y is a scalar operator. Wc conclude that the mapping sending {A,B) S Pp to the 
permutation tta.b G 5*6 is one-to-one. This proves the first assertion. 

The proof of the second assertion is based on matrices U and V from the proof of Theorem [19] and the invariants 
{Jf^,jf), i = 1,2,3, defined in Section [iTj Let us write the product states G kerp, (i = 1,...,6) as Itpi) = 
^ \Xi)- The first five and \xi) are determined directly by U and V, respectively. One has to compute the 
sixth product state in ker p. We have done this in the proof of the theorem mentioned above. Finally, we have verified 
that the 720 (ordered) quintuples that one can construct from the 6 points Pi have generically different values for the 
invariants (J/^, Jf). This completes the proof. □ 

For special states p € £4, the stabilizer may be nontrivial. For instance, the stabilizer for the Pyramid example 
is the alternating group A^ of order 60 which permutes transitively the six product states in the kernel. The original 
definition of this example exhibits only the dihedral group of order 10 as symmetries of a regular pentagonal 
pyramid. A more symmetric realization is given in the recent paper [37j where the full group of symmetries, A5, can 
be realized by local unitary operations. 

For the Tiles example the stabilizer is a group of order 12 isomorphic to the alternating group A4. It also permutes 
transitively all six product states in the kernel of p. Since this stabilizer is finite, it may be conjugated into the 
maximal compact subgroup of PGL which is just the image of the local unitary group. In this way we can obtain a 
more symmetric realization of Tiles. To be concrete, let us consider the following two 3x6 matrices 



U = 



a 1 a 0-1 
1 a -1 a 
1 a -1 a 



V = 



-10 a 1 a 
a -1 a 1 
a -1 a 1 



(121) 



where a = Define the pure product state {k = 1,...,6), as the tensor product of the kth. columns 

of U and V. These states are linearly dependent since the first three and the last three have the same sum. A 
computation shows that the Tiles quintuple IV'fc), k = 0, ...,4, and the quintuple (IV'o)' l''/'4)i IV's) JV'2)) have 
the same invariants. Hence, they are BP-equivalent. The advantage of this new realization of Tiles is that its 
symmetry group (i.e., the stablizcr) is now evident. The symmetry operations arc given by local unitary operations. 
For instance, the multiplication of U and V by diag(— 1, 1, 1) acts on the 6 product states as the permutation (03)(25). 
The multiplication of U by the cyclic matrix Z with first row [0 1] and the simultaneous multiplication of V by 
act as the permutation (012)(345). One can easily construct the unique PPTES p of rank 4 whose kernel is the 
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subspace spanned by the IV"!-)- Up to normalization, p is given by Eq. (I96|) with 
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(122) 



It is now easy to verify that the above two local unitary transformations commute with p. 



VI. PHYSICAL APPLICATIONS 



In this section we demonstrate a few applications of our results in previous sections. 



First, apart from existing numerical tests such as the semidcfinitc programming |ll| . Theorem 1251 analytically 
provides the first bipartite system with specified dimensions and rank, in which all PPTES p can be systematically 
built and characterized. The procedure is as follows. Given a 4-dimensional subspace, one can readily obtain its 
5-dimensional orthogonal complement. Numerically it is possible to find the product states of this subspace. By 
Theorems [TOl and it suffices to consider the case in which this subspace contains exactly 6 product states Xi), 
which are required to be in general position. Next by using Eq. ([2|) and condition (b) of Theorem [191 '^s can compute 
the invariants (Jj*, J^, J^, Jj^, J^, J^) for any quintuple selected from these 6 product states. It is necessary that 
all these invariants be real and that their symbols be of UPB type. Then we compare them with Eq. ([7^ and derive 
the parameters 7a,b,^a,b G (0j''''/2) of a UPB quintuple |ai,6i), see Lemma [T51 Next, we find the ILOs A,B such 
that A ® B\ai, bi) oc \4)i, Xi)- Finally, we build the PPTES A (g) B{I - Yli=o bi){ai,bi\)A'' « B\ which is unique by 
Theorem [Ml This is also succintly stated in Theorem [^ 

On the other hand it is known that the PPT condition is necessary and sufficient for detecting any 2x2 and 2x3 
separable states. Likewise, Theorems [T9l a nd [25l essentially outperform all existing criteria, such as the range criterion 
[27| and the covariance matrix criterion [13|, which can only detect some special 3x3 PPTES of rank 4. It also 
follows easily from Theorem [^ that 



Corollary 28 In a two-qutrit system, the partial transpose of a PPTES of rank 4 has also rank 4- 

Second, we claim that in a two-qutrit system all PPTES of rank 4 are extreme points of the set 5ppT of PPT states. 
(Here our states are assumed to be normalized.) According to the definition in Sec. I, a PPTES p is not extreme 
if and only if it is the midpoint of the segment joining two different PPT states. Let us prove our claim. Suppose 
that 2p = pi + p2 where pi and p2 are two distinct PPT states. Since there is no product state in the range of p, 
the same is true for pi and p2- Thus they are PPTES, and rankpi = rank/92 = 4 since there is no PPTES of rank 2 
or 3 H^. Hence, these three states must have the same range, and so pi ^ p2 = p by Theorem [Mj Thus we have a 
contradiction. Hence p is always both extreme and edge PPTES [l^. The latter statement readily follows from the 
definition of the edge PPTES p, which contains no product state \a, h) G Ti-ip) such that \a* ,h) £ 'R-{p^) [s^ . 

Third, we can systematically build more PPTES and detect them in experiment. By following the technique in 
[3^ , the entanglement witness of the PPTES a of rank 4 has the form 

= P-el, 

e = inf (e,/|P|e,/), (123) 

ie,/) 

where P is the projector onto the kernel of a. The numerical estimation of some states a is also available, as well as 
their experimental realization in (l6j . In this sense, we can detect any two-qutrit PPTES of rank 4 effectively. 

On the other hand it has been proved that any PPTES has the form p = pps + (1 — p)<y where ps is a separable 
state and a is an edge PPTES [Sj]. In the second item, we have shown that any PPTES of rank 4 is an extreme 
and edge PPTES [l9j . So the PPTES p can be characterized when a has rank 4 and the perturbation p is small, 
by using the entanglement witness in Eq. (jl23p . Building the entanglement witness for arbitrary p requires the 
characterization of separable ps, which is still an open problem. Nevertheless, if a PPTES is the convex sum of a few 
PPT states of rank at most 4, it is possibly characterized through the results in this paper. 

Fourth, we claim that no PPTES p g £4 is symmetric, that is, 'R-{p) is not contained in the space spanned by the 
and the + with i > j. To prove this, we use the expression in Eq. (|108p . Assume there is an ILO B = [bij] 
such that I (E) B p I (g) B'' is symmetric. One can readily show that 611 = 621 = = ^22 = 0, which contradicts the 
assumption det B 7^ 0. Hence the simplest symmetric PPTES must have rank at least 5. Such states indeed exist, 
e.g., the state pbea in fssj . 
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Finally, up to ILO, the range of any 3x3 PPTES of rank 4 is orthogonal to a S-dimensional subspacc spanned 
by a UPB. Hence the product vectors of a UPB cannot be distinguished by LOCC 0- In other words, the LOCC- 
indistinguishable nonlocality of the complementary subspace is a deterministic feature of any PPTES of rank 4. 



VII. CONCLUSIONS 



In this paper we have shown that any two-qutrit PPT entangled state of rank 4 is the normalization of 

4 

{A ® B){I - b^){a^M\){A ® B)\ (124) 

i=0 

where A^B are invertible operators and the five product states {\ai,hi)} form an (orthogonal) UPB. Moreover, this 
is the only PPT entangled state among the states having the same range. The 5-dimcnsional subspace spanned by a 
UPB contains exactly 6 product states (up to scalar multiple). We have shown that any 5 of them can be converted 
to a UPB quintuple by a biprojective transformation. The result has been demonstrated on two well-known examples 
of UPB, the Pyramid and Tiles UPB 0]. Therefore we have systematically characterized all PPTES in this system. 

Furthermore we have characterized the separable two-qutrit states of rank 4 whose kernel contains either infinitely 
many product states or exactly 6 product states but not in general position. The next goal in the future is to extend 
our results to higher dimensional PPT states, entangled or separable. 

On the other hand, we have proposed a method of determining the BP-equivalence between two quintuples of 
product states of two qutrits in general position. Apart from the derivation of the main result, the method has also 
been applied to classify the 5-dimcnsional subspaces via their intersection with the set of product states. In particular 
all 11 partitions of the integer 6 occur as the intersection patterns of 'P^C\Ti2.2- These results are useful to characterize 
the distillability of NPT states of rank 4, which is another interesting open problem proposed in Q- 
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